The Varieties of Dynamic Multiscaling in Fluid Turbulence 



cn 
o 
o 

(N 



U 



> 

cn 
o 

On 
O 

m 
o 



> 
X 



Dhrubaditya Mitra and Rahul Pandi10 
Centre for Condensed Matter Theory, Department of Physics, 
Indian Institute of Science, Bangalore 560012, India 

We show that difFerent ways of extracting time scales from time-dependent velocity structure 
functions lead to different dynamic-multiscaling exponents in fluid turbulence. These exponents 
are related to equal-time multiscaling exponents by different classes of bridge relations which we 
derive. We check this explicitly by detailed numerical simulations of the GOY shell model for 
fluid turbulence. Our results can be generalized to any system in which both equal-time and time- 
dependent structure functions show multiscaling. 

PACS numbers: 47.27.i, 47.53.-|-n 



The dynamic scaling of time-dependent correlation 
functions in the vicinity of a critical point was under- 
stood soon after the scaling of equal-time correlations [ij . 
By contrast, the development of an understanding of the 
dynamic mutiscaling of time-dependent velocity struc- 
ture functions in homogeneous, isotropic fluid turbulence 
is still continuing; and studies of it lag far behind their 
analogs for the multiscaling of equal-time velocity struc- 
ture functions 0- There are three major reasons for 
this: (1) The multiscaling of equal-time velocity structure 
functions in fluid turbulence is far more complex than the 
scaling of equal-time correlation functions in critical phe- 
nomena 0- (2) The dynamic scaling of Eulerian- velocity 
structure functions is dominated by sweeping effects that 
relate temporal and spatial scales linearly and thus lead 
to a trivial dynamic-scaling exponent Z£ = \^ where the 
subscript £ stands for Eulerian. (3) Even if this dominant 
temporal scaling because of sweeping effects is removed 
(see below), time-dependent velocity structure functions 
do not have simple scaling forms. As has been recog- 
nized in Ref. Q, in the fluid-turbulence context, an infin- 
ity of dynamic-multiscaling exponents is required. These 
are related to the equal-time multiscaling exponents by 
bridge relations, one class of which were obtained in 
Ref. [3- In the forced-Burgers-turbulence context a few 
bridge relations of another class were obtained in Refs. |3| 
and 0- If the bridge relations of Refs. i and 110 are 
compared naively, then they disagree with each other. 
However, the crucial point about dynamic multiscaling, 
not enunciated clearly hitherto, though partially implicit 
in Refs. is that different ways of extracting time 

scales from time-dependent velocity structure functions 
yield different dynamic-multiscaling exponents that are 
related to the equal-time multiscaling exponents by dif- 
ferent classes of bridge relations. We systematize such 
bridge relations by distinguishing three types of methods 
that can be used to extract time scales; these are based, 
respectively, on integral I, derivative D, and exit-time E 
scales. We then derive the bridge relations for dynamic- 
multiscaling exponents for these three methods. Finally 
we check by an extensive numerical simulation that such 
bridge relations are satisfied in the GOY shell model for 



fluid turbulence. 

To proceed further let us recall that in homogeneous, 
isotropic turbulence, the equal-time, order-p, velocity 
structure function iSp (f) = {[5v\\{x,t,t)Y) ^ for rjd <€. 
^ < L, where Sv\i{x,t,£) = [v{x + It) - v{x,t)] ■ {£/£), 
v{x, t) is the fluid velocity at point x and time t, L is the 
large spatial scale at which energy is injected into the 
system, ry^ is the small length scale at which viscous dis- 
sipation becomes significant, C,p is the order-p, equal-time 
multiscaling exponent, and the angular brackets denote 
an average over the statistical steady state of the tur- 
bulent fluid. The 1941 theory (K41) of Kolmogorov % 
yields the simple scaling result C,!^^^ = p/3. However, ex- 
periments and simulations indicate multiscaling, i.e., (p is 
a nonlinear, convex function of p; and the von-Karman- 
Howarth relation yields ^3 = 1- To study dynamic 
multiscaling we use the longitudinal, time-dependent, 
order-p structure function 

Tp{t,{ti,...,tj,}) = {[5v^\{x,ti,e)...5v\\{x,tp,i)]). (1) 
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restrict ourselves to the simple case ti 
and — tq^2 — ■ ■ ■ — tp = 0, for notational simplic- 
ity write J^p{£,t), and suppress the q dependence which 
should not affect dynamic-multiscaling exponents (see 
below). To remove the sweeping effects mentioned above, 
we must of course use quasi-Lagrangian 0, Q or La- 
grangian |^ velocities in Eq. , but we do not show this 
explicitly here for notational convenience. Given J-p{£, t), 
we can extract a characteristic time scale Tp{£) in sev- 
eral different ways, as we describe later. The dynamic- 
multiscaling ansatz Tp{£) ^ £^p can now be used to de- 
termine the order-p dynamic-multiscaling exponents Zp. 
Furthermore, a naive extension of K41 to dynamic scal- 
ing 13 yields z^4i ^ ^^kh ^ 3/8 for all p. 

In the multifractal model 2] the velocity of a turbulent 
flow is assumed to possess a range of universal scaling 
exponents h ^ X = {hmim hmax) ■ For each h in this 
range, there exists a set S/i C of fractal dimension 
Dih), such that cx (f )'* for r G , with vl the 
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velocity at the forcing scale L, whence 



Sp{i) _ {SvP{i)) 



oc 



d^^{h){- 



\Z{h) 



(2) 



where Z{h) = [ph + 3 — D{h)], the measure d^{h) gives 
the weight of the fractal sets, and a saddle-point eval- 
uation of the integral yields (p = inf[Z(ft,)]. The ph 

h 

term in Z{h) comes from p factors of {£/L) in Eq. ||2J); 
the 3 — D(h) term comes from an additional factor of 
{£/L)^^^^'^\ which is the probability of being within a 
distance £ of the set E/j of dimension D(h) that is 
embedded in three dimensions. Similarly for the time- 
dependent structure function 



£ 



Tp^h 



(3) 



where t?^'''(^r^) has a characteristic decay time Tp ^ ~ 

£/6v{£) - £^-^, and ^^■''(0) = 1. If t'^^'^-^^QP^^dt 
exists, we can define the order-p, degree-M, integral time 
scale 



'^p,M (^) 



Sp{£) 



Tp{£,t)t^^'^~^Ut 



(l/M) 



(4) 



We can now define the integral dynamic-multiscaling ex- 
ponents via M ~ £'^v,M ^ By substituting the mul- 



tifractal form © in Eq. , computing the time integral 
first, and then performing the integration over the mul- 
tifractal measure by the saddle-point method, we obtain 
the integral bridge relations 



4m = 1 + [Cp-M - Q/M, 



(5) 



which was first obtained in Ref. 3]. Likewise, if 
d'^' exists, we can define the order-p, degree-M, 



at 

derivative time scale 



Sp{£) dt^ 



Tp{£,t) 



i-l/M) 



(6) 



the derivative dynamic-multiscaling exponents z via 

Tp^j ~ ^^p *f, and thence obtain the derivative bridge 
relation 



= 1 + [Cp - Cp+m]/M. 



(7) 



Such derivative bridge relations, for the special cases (a) 
p = 2, M = 1 and (b) p = 2, M = 2, were first obtained in 
the forced-Burgers-turbulence context in Refs. H andd, 
respectively [l^, without using quasi-Lagrangian veloc- 
ities but by using other methods to suppress sweeping 
effects. Case (a) yields the interesting result z^^^ = C2, 
since Cs 1- Both relations lO and Q reduce to 
Zp'^^ = 2/3 if we assume K41 scaling for the equal-time 
structure functions. 



If we consider n non-zero time arguments for the 
structure function, Tp^n{£, ^i, • . . , tn, ■ ■ ■ , . . . , 0), which 
we denote by Tp^n{£, ti, ■ ■ ■ , tn) for notational simplicity, 
we can define the integral time scale, TpM^nW = 

VSp{t) Jo ''P\^^ ''1' ■ • ■ ' ''nji'i "■''1 • ■ ■ l-n "'•nj i 

and the derivative time scale, = 
[ Sp(e) dTp' ' ' ' p(^'*l' ■ ■ ■ : *n)|ti=o,...,t„=o]~"^^*-^'^"\ 

where M = From these we can obtain, as 

above, two generalized bridge relations : 



4,M,n = 1 + (Cp-nM - C,p) / (uM); 
Zp,M.n = 1 + (Cp - Cp+um )/ (nM). 



(8) 



In the rest of this paper, we study time-dependent 
structure functions of the GOY shell model for fluid tur- 
bulence HO [13: 



+ iykl)un = i(a„u„+iw„+2 -I- 

6„U„_iU„+i + C„M„„i-U„_2)* + fn- 



(9) 



Here the complex, scalar velocity u„, for the shell n, 
depends on the one-dimensional, logarithmically spaced 
wavevectors fc„ = A:o2", complex conjugation is denoted 
by *, and the coefficients a„ = fc„, 6„ = — (5fc„_i, and 
On = —(1 — (5)fc„_2, with 6 — 1/2, are chosen to con- 
serve the shell-model analogs of energy and helicity in 
the inviscid, unforced limit. By construction, the veloc- 
ity in a given shell is affected directly only by velocities 
in nearest- and next-nearest-neighbor shells. By contrast, 
all Fourier modes of the velocity field interact with each 
other in the Navier-Stokes equation as can be seen eas- 
ily by writing it in wave-vector space. Thus the GOY 
shell model does not have the sweeping effect by which 
modes (eddies) corresponding to the largest length scales 
affect all those at smaller length scales directly. Hence 
it has been suggested that the GOY shell model should 
be thought of as a model for quasi-Lagrangian veloci- 
ties . We might anticipate therefore that GOY- model 
structure functions should not have the trivial dynamic 
scaling associated with Eulerian velocities; we show this 
explicitly below. We integrate the GOY mod el (Et by 
using the slaved, Adams-Bashforth scheme [l^ll5l|. and 
22 shells (1 < n < 22), with /„ = for n > 2 and 
fi — {1 + i) X 5 X 10"'^ (Table nj. The equal-time struc- 
ture function of order-p and the associated exponent is 
defined by 5p(fc„) — (|u„|p) ^ kn'''' ■ However, the static 
solution of Eq. exhibits a 3-cycle with the shell index 
n, which is effectively filtered out 0| if 

(|3(u„+2M„+iMn - (1/4)m„_im„'u„+i)|^'/^) ~ kn'"" , to de- 
termine (^p. These exponents are in close agreement with 
those found for hornogeneous, isotropic fluid turbulence 
in three dimension |l2l |. Data for the exponents (p from 
our calculations are given in Table We analyse the 
velocity (u„(t)) time-series for n = 4 to 13, which cor- 
responds to wave-vectors well within the inertial range. 



3 



V 


it 


A 


Urms 


Rex 
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Tav 


10-^ 


2 X 10"'' 


0.7 


0.35 


2 X lO'' 


6.3 




5 X W^St 


IOVl 



TABLE I: Viscosity v, the time-step St, Taylor microscale 
^ = (I]„ l^"l^/^n/I]„ the root-mean-square 

velocity Urms = [2 |itn|^/fcn/(27rfci)]"'"''^, the Taylor- 
microscale Reynolds number Rex = \urms/v, the integral 
scale Lint = hn|^/fcJ;)/(^„ |M„|^/fc„), and the box-size 
eddy turnover time = Lint/Urms, that we use in our nu- 
merical simulation of the GOY shell model. Data from the 
first Ttr time steps are discarded so that transients can die 
down. We then average our data for time-dependent struc- 
ture functions for an averaging time Tav 



The smaller the wave- vector fc„ the slower is the evolu- 
tion of Unit), so it is important to use different temporal 
sampling rates for velocities in different shells. We use 
sampling rates of 50 x 6t for 4 < n < 8 and 10 x St for 
9 < ri < 13, respectively. 

For the GOY shell model we use the normalized, 
order-p, complex, time-dependent structure function, 
fp{n,t) = {[un{0)u*„{t)]P/'^)/Sp{kn) , which has both real 
and imaginary parts. The representative plot of Fig. ^ 
shows that the imaginary part of fp{n,t) is negligibly 
small compared to its real part. Hence we work with the 
real part of fp{n,t), i.e., Fp{n,t) = 5ft[/p(n,t)]. 

Integral and derivative time scales can be defined for 
the shell model Q as in Eqs. Q and ©. We now 
concentrate on the integral time scale with M = 1, 
Tp i{n, tu) = /q" Fp{n, t)dt, the derivative time scale with 

M = 2, Tp^2 = [^-|F^I*=o]"^^^ and the associated 
dynamic-multiscaling exponents defined via ^{n, t^) ^ 

-z' 

kn '''^ and Tp2{'^) ~ ^''^ ■ In principle we should use 
t„ — > cxD but, since it is not possible to obtain Fp{n,t) 
accurately for large t, we select an upper cut-off such 
that Fp{n,tu) = a, where, for all n and p, we choose 
a = 0.7 in the results we report. We have checked that 
our results do not change if we use 0.3 < a < 0.8. The 
slope of a log-log plot of Tp i{n) versus fc„ now yields 
Zp I (Fig. Hand Table |Tll). Preliminary data for ^ were 
reported by us in Rcf. 9J. 

For extracting the derivative scale extend 
Fp{n,t) to negative t via Fp{n, —t) = Fp{n,t) and use 
a centered, sixth-order, finite-difference scheme to find 



j^Fp{n,t)\t=Q. A log-log plot of T,^ 



p2('^) versus fc„ now 



yields the exponent (Fig. ^ and Table llT|) . 

In Ref. 13] dynamic-multiscaling exponents were ex- 
tracted not from time-dependent structure functions but 
by using the following exit-time algorithm: We define the 
decorrelation time for shell n, at time ti, to be Ti{n), such 
that, \un{ti)\\un{U + Ti)\ > A±i|u„(t,)P,withO< A< 1. 
The exit-time scale of order-p and degree-M for the shell 



rpE 

^p,M 



lim 

N^oc 



N 



N 



(10) 



where the last proportionality follows from the dynamic- 
multiscaling ansatz. In practice we cannot of course take 
the limit A'^ — > 00; in a typical run of length Tav (Ta- 
ble A^ ~ 10^. By suitably adapting the multifractal 
formalism used above, we get the exit-time bridge rela- 



tion Zpj^j = 1-1- [Cp-M — Cp]/-^j obtained in Ref. 
only for M = 1. Dynamic-multiscaling exponents ob- 
tained via this exit-time algorithm are shown for M = 1 
and AI = —2 in Table HTl The exit-time bridge relations 
for M > are the analogs of the integral-time bridge 
relation (|SJ and those for M < are the analogs of the 
derivative-time bridge relation {T)) . We have checked that 
our results do not depend on A for 0.3 < A < 0.8. 

Our numerical results for the equal-time exponents Cp 
(Column 2), the integral-time exponents Zp i (Columns 
3 and 4), the derivative-time exponents Zp2 (Columns 
6 and 7), and the exit-time exponents Zp^i and Zp _2 
(Columns 5 and 8, respectively) for 1 < p < 6 are given 
in Tabled] The agreement of the exponents in Columns 
3 and 4 shows that the bridge relation Q is satisfied 
(within error bars). Likewise, a comparison of Columns 
6 and 7 shows that the bridge relation O is satisfied. By 
comparing Columns 4 and 5 we see that the integral-time 
exponent Zp i is the same as the exit-time exponent Zpi ; 
similarly. Columns 7 and 8 show that the derivative-time 
exponent Zp2 is the same as the exit-time exponent Zp_2- 
The relation Z21 = C2 mentioned above |^ is not mean- 
ingful in the GOY model since dFp{n, t)/dt\t=o vanishes, 
at least at the level of accuracy of our numerical study. 

We have obtained 50 different values of each of the 
dynamic-multiscaling exponents from 50 different initial 
conditions. For each of these initial conditions time- 
averaging is done over a time Tat, (Table P) which is 
larger than the averaging time of Ref. ^3 a factor 
of about 10^. The means of these 50 values for each 
of the dynamic-multiscaling exponents are shown in Ta- 
ble^ and the standard deviation yields the error. This 
averaging is another way of removing the effects of the 
3-cycle mentioned above. 

We have shown systematically how different ways of ex- 
tracting time scales from time-dependent velocity struc- 
ture functions or time series can lead to different sets 
of dynamic-multiscaling exponents, which are related in 
turn to the equal-time multiscaling exponents ^p by dif- 
ferent classes of bridge relations. Our extensive numerical 
study of the GOY shell model for fluid turbulence verifies 
explicitly that such bridge relations hold. Experimental 
studies of Lagrangian quantities in turbulence have been 
increasing over the past few years We hope our 

work will stimulate studies of dynamic multiscaling in 
such experiments. Furthermore, the sorts of bridge rela- 
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FIG. 1: Plots of real (ai) and imaginary (a2) parts of the time-dependent structure function fp{n,t) for the GOY shell model 
for shell number n = 9 and order p — 5 versus time t/rL, where tl is the box-size eddy turnover time (Table HJ. Note that 
Q[fp{n,t)] is negligibly small compared to Fp(n,t) = lR\fp(n,t)]. (b) Fp(n,t) versus t/rh for p = 5 and n = 5, 7, and 9. 
Representative log-log plots (base 10) of the integral (c) and derivative (d) time scales T/_i(n) and T^2{'n) versus kn\ the slopes 
of the linear least-square fits in (c) and (d) yield the dynamic exponents i and 23^2, respectively. 



order (p) 


Cp 


4i[Eq-0] 








Zp,2 


^p,-2 


1 


0.3777 ± 0.0001 


0.6221 ± 0.0001 


0.60 ± 0.02 


0.603 ± 0.007 


0.6820 ± 0.0001 


0.70 ± 0.02 


0.677 ±0.001 


2 


0.7091 ± 0.0001 


0.6686 ± 0.0002 


0.67 ± 0.02 


0.661± 0.007 


0.7081 ± 0.0002 


0.71 ± 0.01 


0.719 ±0.004 


3 


1.0059 ± 0.0001 


0.7030 ± 0.0002 


0.701 ± 0.009 


0.708 ± 0.001 


0.7310 ± 0.0002 


0.73 ± 0.01 


0.739 ±0.006 


4 


1.2762 ± 0.0002 


0.7298 ± 0.0003 


0.727 ± 0.007 


0.74±0.01 


0.7509 ± 0.0003 


0.744 ± 0.009 


0.758 ±0.006 


5 


1.5254 ± 0.0005 


0.7511 ± 0.0007 


0.759 ± 0.009 


0.77± 0.01 


0.7684 ± 0.0007 


0.756 ± 0.009 


0.778 ±0.003 


6 


1.757 ± 0.001 


0.768 ± 0.002 


0.77 ± 0.01 


0.79± 0.01 


0.7836 ± 0.002 


0.764 ± 0.009 


0.797 ±0.0008 



TABLE II; Order— p (Column 1) multiscaling exponents for 1 < p < 6 from our simulations of the GOY model; equal-time 
exponents (^p (Column 2), integral-scale dynamic-multiscaling exponent of degree-1 (Column 3) from the bridge relation @ 
and the values of (^p in Column 1, Zp i from our calculation using time-dependent structure functions (Column 4), the exit-time 
exponents of order 1 Zp i (Column 5), the derivative-time exponents z^2 (Column 6) from the bridge relation Q and the values 
of d^p in Column 1, 2^2 from our calculation using time-dependent structure function (Column 7) and the exit-time exponent 
of order -2, Zp_2 (Column 8). The error estimates are obtained as described in the text. 



tions we have discussed here must also hold in other prob- 
lems with multiscaling of equal-time and time-dependent 
structure functions or correlation functions. Passive- 
scalar and magnetohydrodynamic turbulence are two ob- 
vious examples which we will report on elsewhere |l7j . 
Numerical studies of time-dependent, quasi-Lagrangian- 
velocity structure functions in the Navier-Stokes equa- 
tion, already under way, will also be discussed elsewhere. 
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